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Abstract. We consider bivariate polynomials orthogonal on the bicircle with 
respect to a positive linear functional. The lexicographical and reverse lexico- 
graphical orderings are used to order the monomials. Recurrence formulas are 
derived between the polynomials of different degrees. These formulas link the 
orthogonal polynomials constructed using the lexicographical ordering with 
those constructed using the reverse lexicographical ordering. Relations be- 
tween the coefficients in the recurrence formulas are derived and used to give 
necessary and sufficient conditions for the existence of a positive linear func- 
tional. These results are then used to construct a class of two variable measures 
supported on the bicircle that are given by one over the magnitude squared of 
a stable polynomial. Applications to Fejer-Ricsz factorization arc also given. 
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1. Introduction 

Bivariate polynomials orthogonal on the bicircle have been investigated mostly 
in the electrical engineering community in relation to the design of stable recursive 
filters for two-dimensional filtering. In particular we note the work of Genin and 
Kamp [7J who were interested in the following problem. Given any two variable 
polynomial q(z, w) with g(0,0) ^ 0, let ak,i(z,w) be its planar least squares inverse 
polynomial of degree (fc, I); i.e., a^i minimizes the mean quadratic value of 1 — a& \q 
on the bicircle. What properties does a^j have? At the time it was conjectured 
the minimizing polynomials were stable, i.e. ak.i(z,w) ^ 0, \z\ < 1, \w\ < 1, 
which they showed was false. Their investigation was carried further by Delsarte, 
Genin and Kamp [4] who developed the connection between these polynomials and 
matrix polynomials orthogonal on the unit circle [3]. In the development of this 
connection these authors were lead to examine moment matrices that were block 
Toeplitz matrices where each block entry is itself a Toeplitz matrix. Such structured 
matrices are called doubly Toeplitz matrices and arise naturally in the bivariate 
trigonometric moment problem. These types of matrices arose more recently in 
the work of Geronimo and Woerdeman [8] in their investigation of the bivariate 
Fejer-Riesz factorization theorem. These authors were able to resolve the question 
when a strictly positive bivariate trigonometric polynomial of a certain degree can 
be written as the magnitude squared of a stable polynomial of the same degree. 
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In this work the authors used the fact that the theory of orthogonal polynomials 
on the unit circle provides a proof of the one variable Fejer-Riesz theorem which 
does not use the fundamental theorem of algebra. We intend here to continue to 
investigate the properties of bivariate polynomials orthogonal on the bicircle and 
clarify their role in the Fejer-Riesz Theorem. 

A major difficulty encountered in the theory of orthogonal polynomials of more 
than one variable is which monomial ordering to use. For bivariate real orthogonal 
polynomials the preferred ordering is the total degree ordering which is the one set 
by Jackson [T3] . For polynomials with the same total degree the ordering is lexico- 
graphical. As noted in Dclgado et al [2] in their study of orthogonal polynomials 
associated with doubly Hankel matrices there is a good reason for choosing this 
ordering which is that if new orthogonal polynomials of higher degree are to be 
constructed then their orthogonality relations will not affect the relations govern- 
ing the lower degree polynomials. However in order for the moment matrix to be 
doubly Toeplitz the monomial orderings that need to be used are lexicographical 
and reverse lexicographical. 

We begin in Section 2 by considering finite dimensional subspaces spanned by the 
monomials z l w J , |t| < n, \j\ < m, and exhibiting the connection between positive 
linear functionals defined on this space and positive definite doubly Toeplitz matri- 
ces. We then introduce certain matrix orthogonal polynomials and show how they 
give the Cholesky factors for the inverse of the doubly Toeplitz matrices considered 
above. The results in [8] show that these polynomials play a role in the parametric 
moment problem. In Section 3 we construct two variable orthogonal polynomials, 
where the monomials are ordered according to the lexicographical ordering. When 
these polynomials are organized into vector orthogonal polynomials they can be 
related to the matrix orthogonal polynomials constructed previously. From this 
relation it is shown that these vector polynomials are the minimizers of a certain 
quadratic functional. Using the orthogonality relation, recurrence relations satisfied 
by the vector polynomials and their counterparts in the reverse lexicographical or- 
dering are derived and relations between these recurrence coefficients are exhibited. 
In Section 4 a number of Christoffcl-Darboux like formulas are derived. In Section 
5 we use the relations between the coefficients derived in Section 3 to develop an 
algorithm to construct the coefficients in the recurrence formulas at a particular 
level (n, m) say, in terms of the coefficients at the previous levels plus a certain 
number of unknowns. The collection of these unknowns is in one to one correspon- 
dence with the number of moments needed to construct the vector polynomials up 
to level (n, m). This is used in Section 6 to construct a positive linear functional 
from the recurrence coefficients. The construction allows us to find necessary and 
sufficient conditions on the recurrence coefficients for the existence of a positive 
linear functional which is in one to one correspondence with the set of positive defi- 
nite doubly Toeplitz matrices. In Section 7 wc examine conditions under which the 
linear functional can be represented as a positive measure supported on the bicircle 
having the form of one over the magnitude squared of a stable polynomial. This 
gives a new proof of the Fejer-Riesz result of [8]. Finally in Section 8 examples are 
given that illustrate various aspects of the theory developed. 
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2. Positive linear functionals and Doubly Toeplitz matrices 



In this section we consider moment matrices associated with the lexicographical 
ordering which is defined by 

(k,£) < lox (h,£i) k < fci or (k = h and £ < £{), 

and the reverse lexicographical ordering defined by 

{k,£) < rov iox {kuh) (£,k) < lox (£i,fci). 

Both of these orderings are linear orders and in addition they satisfy 

(k, £) < (m, n) (k + p, £ + q) < (m + p,n + q). 

In such a case, one may associate a halfspace with the ordering which is defined by 
{(k, I) : (0, 0) < (k, I)}. In the case of the lexicographical ordering we shall denote 
the associated halfspace by H and refer to it as the standard halfspace. In the case 
of the reverse lexicographical ordering we shall denote the associated halfspace by 
H. Instead of starting with the ordering, one may also start with a halfspace H of 
Z 2 (i.e., a set H satisfying H + H <zH,HC\ (-H) = 0, # U (-H) U {(0, 0)} = I?) 
and define an ordering via 

(k, I) < k {k u h) (fci -k,h-l)€ H. 



We shall refer to the order <^ as the order associated with H 



Note that the 

lexicographical and reverse lexicographical orderings do not respect total degree. 

Let J^™'" 1 denote the bivariate Laurent linear subspace spanjz 2 ^ , — n < i < 
n, — m < j < to}. Let C n ,m be a linear functional defined on rj"' m by 



W 



We will call c,j- the (i, j) moment of £„ lTO and j£ n ,m a moment functional. If we 
form the (n+ l)(m+ 1) x (n+ 1)(to+ 1) matrix C„ iTO for C n ,m in the lexicographical 
ordering then, as noted in the introduction, it has the special block Toeplitz form 

Cq C-i ■ ■ ■ C- n 



C\ Co 



C_ 



n+l 



(2.1) 



C n C n _i • • • Co 
where each Q is an (to + 1) x (to + 1) Toeplitz matrix as follows 



Ci.. 



Cifi 



i = —n. 



,71. 



(2.2) 



Thus C„ jTra has a doubly Toeplitz structure. If the reverse lexicographical ordering 
is used in place of the lexicographical ordering we obtain another moment matrix 
C n , m where the roles of n and m are interchanged. 

Let us introduce the notion of centro-transpose symmetry. We denote the trans- 
pose of a matrix A by A T . A square matrix A is said to be centro-transpose 
symmetric if JAJ = A T where J is the matrix with l's on the antidiagonal and 
zeros elsewhere. Note that a Toeplitz matrix is centro-transpose symmetric. We 
have the following useful lemmas which characterize Toeplitz and doubly Toeplitz 
matrices in terms of centro-transpose symmetry. 
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Lemma 2.1. An (n + 1) X (n+ 1) matrix A — (cii,j)fj—o is Toeplitz if and only if 
both A and A := (di^fj—o are centro-transpose symmetric. 

Proof. Notice that JAJ = A T is equivalent to a n ^i_ n ^j — a^j, < i,j < n. 
Similarly, the centro-transpose symmetry of A is equivalent to a n _i_j n _i_j = atj j, 
< i,j < n— 1. But then 

Oi+lJ+l = On-j-l.n-i-l = < Z, j < n - 1, 

and thus it follows that ^4 is Toeplitz. 

As ^4 and A are Toeplitz, the converse is immediate. □ 

Lemma 2.2. Let A = (Aij), i,j = 1, . . . , k, where each Ai t j is a complex m x to 
matrix. Then A is a doubly Toeplitz matrix if and only if A T = JAJ, Af = J\ A\ J\ 
and A2 = J1A2J1. Here A\ is obtained from A by deleting the last block row and 
column and A2 is obtained from A by removing the last row and column of each 
Aij . The matrices J and J\ are square matrices of appropriate size with ones on 
the antidiagonal and zeros everywhere else. 

Proof. Again the necessary conditions follow from the structure of A. To see the 
converse note that A T = JAJ implies that A 1 ^ i — Ji-^k-i,h-jJi, where J2 is the 
to x to matrix with ones on the reverse diagonal and zeros everywhere else . This 
coupled with the condition on A\ implies that A is a block Toeplitz matrix from 
Lemma [2.11 and J2^4j.jJ2 = Afj. These relations plus the condition on A2 and 
Lemma |2~T1 gives the result. □ 

Remark 2.3. The conclusion of the above lemmas hold if we replace deleting the 
last (last block) row and column by deleting the first (first block) row and column. 

We say that the moment functional £ n , m ■ Yl ' — > C is positive definite or 
positive semidcfinite if 

(\p\ ) > or C n!m (\p\ ) > (2.3) 

for every every nonzero polynomial p £ J}™'™. It follows from a simple quadratic 
form argument that C n ,m is positive definite or positive semidcfinite if and only if 
its moment matrix C n>m is positive definite or positive semidefinite respectively. 
We will say that C is positive definite or positive semidefinite if 

£(H 2 )>0 or £(M 2 )>0 

respectively for all nonzero polynomials. Again these conditions are equivalent to 
the moment matrices C n ^ m being positive definite or positive semidefinite for all 
positive integers n and m. The above discussion leads to, 

Lemma 2.4. Let C„, m be a positive (positive semi-) definite (n + l)(m + 1) x (n + 
l)(m+ 1) matrix given by (|2.ip and (|2.2|) . Then there is a positive (positive semi-) 
definite moment functional : J^ n,m — > C associated with C n , m given by 

Cij = £ n ^ n (z~ l w~ j ) = £ n ,m(z*w J '), -« < i < ti, —m<j<m. 
The converse also holds. 

Let Ilm+i ^ e the set of all (to+ 1) x (m + 1) complex valued matrix polynomials 
of degree n or less, Jlm+i the set of all (to + 1) x (m + 1) complex valued matrix 
polynomials, and M m ' n is the space of to x n matrices. For a matrix M we let M* 
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denote the conjugate transpose (or the adjoint) of M. For a polynomial Q(z, w) we 
let Q^(z, w) denote the polynomial in z^ 1 and w" 1 defined by Q^(z, w) — Q( = , =Y ■ 
If the positive moment functional C n . m : Y[ n ' m — > C is extended to two variable 
polynomials with matrix coefficients in the obvious way, we can associate to it a 
positive matrix function C m : Ilm+i x Ilm+i ^ M m+1 ' m+1 defined by, 

[C m (P(z),Q(z))] ij = C n , m ([P(z,w) Q<<(z,w)].J, l<i,j<m+l (2.4) 

where 



P(z,w) = P(z) 



w" 



and Q(z,w) — Q(z) 



Equation (|2.4[) shows that if C n ,m can be represented as a positive measure \x 
supported on the bicircle then for /anm+lxm + 1 matrix function continuous 
on the unit circle, 



Crn(f) = / f(6)dM m (6), 

J — 7T 

where M m is the m + 1 X m + 1 matrix measure given by, 



dM m (0) 



which shows that M m is Toeplitz. 

Because of the structure of C n ^ m we can associate to C m matrix valued orthogonal 
polynomials in the following manner [3J, [3], [5]. Let {i?™( 2: )}r=o anc ^ {^T( z )}i=o 
be (to + 1) x (to + 1) complex valued matrix polynomials given by 



" 1 " 




' i " 








w ra 







and 

satisfying 
and 



Lr(z)=L™z* + L™„ 1 z*- 1 + -- - , 



i = 0, . . . , n, 

i = 0, . . . , n, 



(2.5) 
(2.6) 
(2.7) 
(2.8) 



respectively, where I m +\ denotes the (to + 1) x (m+ 1) identity matrix. The above 
relations uniquely determine the sequences {i?™}™ =0 an< ^ {^rlr=o U P to a unitary 
factor and this factor will be fixed by requiring i?™ and L™ to be upper triangular 
matrices with positive diagonal entries. We write 

TIT" 
4/ -I ft 



L™(z) = [()••• i™ iy_i •• L 



i,0i 



777+1 
77—1 J 

Z Jm+1 



m+1 



(2.9) 



and 







Z J-m+1 






77—1 T 

Z Jm+1 


= L 




. Lh n (z) _ 




-^777+1 



(2.10) 
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where 



r m T m 

n,n n,n— 1 

n t m 

u -^n-l.n-l 





■ J n,0 



L< 



J 0,0 



(2.11) 



In an analogous fashion write 



Rtf(z} 
R?(z) 



m+1 



zU Im+l] R, 



(2.12) 



where 



R 



-"-0,0 "l.O 

iZf x 











Dm 
pm 



(2.13) 



By lower (respectively upper) Cholesky factor A (respectively £?) of a positive 
definite matrix M we mean 



M = AA^ = BBK 



(2.14) 



where A is a lower triangular matrix with positive diagonal elements and B is an 
upper triangular matrix with positive diagonal elements. With the above we have 
the following well known lemma 15J, 



Lemma 2.5. Let C n ^ m be a positive definite block Toeplitz matrix given by (12.1 
then U is the lower Cholesky factor and R is the upper Cholesky factor of C~ x m . 



Proof. To obtain ([2~T4]) note that ([2T8]> implies that 





( 


Z J-m+l 




Z J-m+l 


\ 






Z Jm+1 




11— 1 T 

Z J-m+l 


C m ((L™),(LZ))=L£ m 
















Im+l 




Im+l 


J 



where I is the (n + l)(m + 1) x (n + l)(m + 1) identity matrix. S ince C n m is 
invertible we find, 

C„ = L L. 



The result for R follows in an analogous manner 
From this formula and (|2. 1 1[) we find, 
LZ(z)=[(L™J)-\0,0,...0 

and 



□ 



Cn,m[ z ™^m+li z ™ , ■ ■ ■ j Im+l] T , (2-15) 



i?™(z) = [I m+1 , zl m+1 , . . . , z"/ m+1 ] C- m [0, 0, ... , 0, (iCJ" 1 ] . (2.16) 

Note that L™ n ' is the lower Cholesky factor of [I m +i , 0, • • • , 0]C~ m [/ m +i , 0, • • • , 0] T 
while Rnn 1S tne upper Cholesky factor of [0, • • • , I m+ i]C~ m [0, • • • , I m +i] T ■ 
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The theory of matrix orthogonal polynomials ([3], [15], [II], [19]) can be applied 
to obtain the recurrence formulas 



A i+1 , m L™ 1 (z)=zL7 l (z)-E i+1 , rn R i (z) * = 0,...,n - 1, 

A < — m 

R%. l (z)A i+ltm =zR?(z)- L i (z)E l+1 . m t = 0,...,n-l, 

where 



(2.17) 



£i+i,m = £ m (zL™, %) = ^™(ir f , {zR?) ] ) (2-18) 



and 

4 — C { vT m T m \ — T m I T m \ ~ 1 

A i+1 . m = C m (R™ +l , (zi?™) f ) = (^r+i.i+i)^ 1 ^™- 



(2.19) 



For a matrix polynomial i? of degree n in z, B(z) = z"Bt(l/z). By multiplying 
the first equation in (|2.17|) on the left by zL™(z)^ and the second equation on the 
right by z_R™(z)^ then integrating we see that 

Ai+i, m A\ +lrn = I rn +\ — Ei + i im Ej +lm , 

The above equations and the properties of Ai + i yTn and Ai + i^ m show that is a 

strictly contractive matrix and that Ai + \ yTn is the upper Cholesky factor of I m +i — 
E i+1 ^ m Ej +l m . Similarly A\ m is the lower Cholesky factor of I m+1 - Ej +l m E i+1 ^ n . 
Furthermore f2~T!)j) and (j^U|) show that 

i+l 

det((L™ lti+1 ?L™ +1 . +1 )- 1 = det(Co) [] det(7 m+1 - E 3 , m E] j. (2.21) 

i=i 

The recurrence formulas (|2 . 1 T[) can be inverted in the following manner. Multiply 
the reverse of the second equation in (|2.17|) on the right by Ei + i_ m to obtain 

E i+ i tm A\ +l m R™ +1 {z) = E i+ i t7n R™(z) - zE l+1 , m E} +1 m L™{z). 

Add this equation to the first equation in (|2.17[) then use (|2.20[) to eliminate A i+ i m 
and A\ , j m to find 

(4 + i, m ) _1 ^iW + E i+1>m (A i+1 , m )- l R? +1 (z) = zLT(z). (2.22) 

In a similar manner we find 

R? +1 {z){A\ +ljm )- 1 + ^W^)- 1 ^!,™ = zJC(«)- (2-23) 

From the recurrence formulas it is not difficult to derive the Christoffcl-Darboux 
formulas [3], 

Rf(z) RT( Zl ) - zz 1 L^( z yLf(z 1 ) = (1 - zzx) ]T LrW^T^i) 

i=0 (2.24) 

L^L^/ - ^i^i)^)* = (1 - zz^RTiz^RTi^ 



These formulas give rise to the matrix Gohbcrg-Semencul formulas |llj.|15j when 
the linear equations obtained by equating like powers of z l z J are put in matrix form. 
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Some properties that follow from the above formulas (3j [Theorems 9, 14 and 15] 
are that R™(z) and L™{z) have empty kernels for \z\ < 1, i.e., 

det(i?™(z)) ? det(L%(z)), \z\ < 1. (2.25) 
Such polynomials are called stable matrix polynomials and if we write 



Wk(z) 



(2.26) 



and 
then 

Furthermore 



C) = ^£ e- ije W k (e ie )d6 



Ci = Cj, \j\<k. (2.27) 



w k = 



(2.28) 



If L™(z) (RT(z)) satisfies (I2.25|) and (|2.27p we will say it is stable and has spectral 
matching (up to level k). Another useful result shown in [3] is 



logdet((L£ M+1 )tL» M+1 )-i = _L £ i og det 



W k (9)d9. (2.29) 



From the stability of Rf +1 and i™ x , (|2~22]l and (|2~23|) give the following formulas 
for the recurrence coefficients Ei + i^ rn : 

E i+1 , m = -(4+l, TO )" li i+l(Q)%i(0)" 1 ^+l,r» (2-30) 

= -A i+ltm LT +1 (0)- 1 RT +1 (0)(A\ +lt J- 1 . 

We also note that L'T'(z) and RT(z) are minimizers of certain quadratic functions. 
To see this denote the set of (to + 1) x (m + 1) hermitian matrices as Herm(m + 1) 
and let M : [] m ,j — > Herm(m + 1) be given by, 

Af[X(z)] - £ m (X,X) - (X(0) + X(0)t), (2.31) 
then Delsarte, Genin and Kamp have shown [3] that for a given degree k, M. 
is minimized by *L%(z)L% m with value (L% m )iL% m . Likewise : U m +i -» 
Herm(m + 1) given by, 

= A»(Xt, Jft) - (X(0) + *(())*), (2.32) 
is minimized by R™ m R™(z) and takes the value R™ m {R™ m y ■ Thus we find, 

{L™ m ) L™ m > (ifc+l, m ) -^fc+l, m (2.33) 

and 

^k,m\^k,m) — -^fc+l.mC^fc+l.m) • (2.34) 

Here A > B for two (m + 1) x (to + 1) matrices means that A — B is positive 
semidefinite. The above discussion leads to Burg's entropy Theorem. Consider the 
class of M m of (to + 1) x (m + 1) matrix Borel measures on the unit circle and 
for each such measure fj, write the Lebesgue decomposition of /i = fi ac + fj, g where 
dn ac /d0 = W{6). Let S™ be the subset of M m such that each /i£S„ m has the 
same Fourier coefficients C i; \i\ < n and £(p) = lndet(W)e£0 > — oo. Then 
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there is a unique measure which maximizes the above entropy function £(/j,) and 
this measure is given by d/j, = W(9)d9 with W{6) = Q™^) -1 where Q™(0) is a 
positive (to + 1) x (m + 1) matrix trigonometric polynomial of degree n. 

This leads to a simple proof of the Matrix Fejer-Reisz factorization Theorem 
(Helson |13j . Dritschel [5], McLean- Woerdeman [16] , Geronimo-Lai [10]) which will 
be useful later. 

Lemma 2.6. Let Q™(#) be a strictly positive (m + 1) x (to+1) matrix trigonometric 
polynomial then Q™{0) = L™(z)(L™(z))\ z = e l9 where L™ is a stable (to + 
1) x (m + 1) matrix polynomial of degree n. Furthermore L™ is given by (|2.15[) . 

Proof. Since Q™(9) is strictly positive we can compute the moments 
Cj = y~ J e~ i: > Q" L (O)" 1 d9 . If we compute the matrix orthogonal polynomials 
associated with these Fourier coefficients we find that W n has spectral matching up 
to n. That is its Fourier coefficients match C for \i\ < n. The Maximum Entropy 
Theorem implies that Q™(0) — W~ x , which gives the result. □ 

The matrix Fejer-Riesz Theorem now follows. 

Theorem 2.7. Let Q™(#) > be a positive (to + 1) x (to+ 1) matrix trigonometric 
polynomial then Q™(0) — P™{z)(P™(z))\ z — e where P™ is an outer (nonzero 
for \z\ < 1) (to + 1) x (to + 1) matrix polynomial. 

Proof. Let Q™ e = el + Q™ e > then Q™ e satisfies the hypotheses of the above 
lemma. Thus Q™ e = P™ e (P™ e ) T . The proof now follows by taking the limit as e 
tends to zero. □ 

It was observed by Delsarte et. al. [3] that if the C& in C n ^ m are centro-transpose 
symmetric then 

{L$L?{z)) T = J m R?{z)Rf J m i = 0,...,n (2.35) 

where J m is the (to+1) x (m + 1) matrix with ones on the reverse diagonal and zeros 
everywhere else. This can easily be seen from formulas (|2 . 15f> and (|2.16|) since in 
this case from Lemma l2"?2"l Cj m = JC n , m J with J the (n+l)(TO+l) x (m + l)(n+l) 
matrix with ones down the anti-diagonal and zeros everywhere else . This leads to 
the following characterization of positive definite doubly Toeplitz matrices in terms 
of certain recurrence coefficients. We will denote by C™ the to x to matrix obtained 
from Co by eliminating the first row and first column of Co. 

Theorem 2.8. Suppose C n ^ rn is positive definite. Then the Fourier coefficients 
Ci, \i\ < n, are centro-transpose symmetric if and only if Ek,m, k = 1, . . . , n, and 
Co are centro-transpose symmetric. Consequently, C n ,m is doubly Toeplitz if and 
only if Ek^i, k = 1, . . . , n, i = m— 1, m, Co and C™ are centro-transpose symmetric. 

Proof. Examining the leading coefficients in (|2.35[) and using the fact that 
and i?™ are upper triangular we find that (see also [4]), (L™) T = J m R"\Jm for 
i — 0, . . . , n. Thus 

L™(zf = J m R™{z)J m , i = 0,...,n (2.36) 
From (|2.18p we find, since J m J m = I m +i, 
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To show the converse note that if -Ei, m is centro-transpose symmetric then from 
(f2T20j) we obtain 

so that 

J ra^i^rnJ va ~ ^i^ra' (2.37) 

Since Co is centro-transpose symmetric and Lo,m(z) is the upper Cholesky factor 
of Co, we see that JmLffJm = R™ 7 . Thus by induction using (|2.17[) we find that 
J m L™(z)J m = R™(z) T . The first part of the result now follows from the spectral 
matching of W n ( 1)2.270 ) and (|2.28|) . The second part of the Theorem follows by 
applying the above argument to C„ jm _i and C™ then using Lemma |2. II . □ 

In the next two sections we present recurrence formulas and an algorithm that 
computes recurrence coefficients for a positive definite doubly Toeplitz matrix. 



3. BlVARIATE ORTHOGONAL POLYNOMIALS 



In this section we examine the properties of two variable orthogonal polynomials 
where the monomial ordering is either lexicographical or reverse lexicographical. 
The study of orthogonal polynomials on the bicircle with this ordering was begun 
by Delsarte et.al. [4] and extended in [8]. Given a positive definite linear functional 
£-n,m '■ II — >■ C we perform the Gram-Schmidt procedure using the lexico- 
graphical ordering and define the orthonormal polynomials tfi l n m (z , w) , < n < 
N, < to < M, < Z < to, by the equations, 



-N.M\ 



A* I 



II .III 



r ) = o, 

) f ) = l, 



< i < n and < j < m or i = n and < j < I, 

(3.1) 



and 



<f>>(z ) w) = k%l n ,z n w l + 



II. I 



(i,j)<lc X (nA) 



(3.2) 



With the convention fc™' , > 0, the above equations uniquely specify 



Poly- 

nomials orthonormal with respect to Cm,m but using the reverse lexicographical 
ordering will be denoted by 4> nm - They are uniquely determined by the above 
relations with the roles of n and m interchanged. 
Set 



where the (m + 1) x (n + 1)(to + 1) matrix K n ^ rn is given by 



^71,171 




- z n w m - 


Am- 1 

n,m 


— K 


z n w m-l 






1 



(3.3) 



K — 



n.m.rn 



h,n,m — 1 

i n,m— 1 
n.m.m— 1 



, n : 

"Vi,,m,0 



1 n—l.m 



1.0,0 

n,m,m— 1 



i 0,0 

C n,m,0 



(3.4) 
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As indicated above denote 



n,m 
in- 1 
n.rn 



W"'Z" 
,..m ~n— 1 



1 



(3.5) 



where the (n+ 1) x (n+ l)(m+ 1) matrix K n , m is given similarly to (|3.4|) with the 
roles of 71 and m interchanged. For the bivariate polynomials <j) l n m (z, w) above we 

define the reverse polynomials <\> l nm {z, w) by the relation 



0' (z,w)=z n w m cl><(l/z,l/w) 



(3.6) 



With this definition <j> n m (z, w) is again a polynomial in z and w, and furthermore 



$„, m (z, W ) ^zW^.Jl/z, !/«>) = 



n,m 

m— 1 

n,m 



n.m 



(3.7) 



An analogous procedure is used to define (f> l nm - 

In order to ease the notation to find recurrence formulas for the vector polyno- 
mials <& n ,m we introduce the inner product, 



(X,Y) =C NM {XY*). 



(3.8) 



Let Yl be the linear span of z l w° , < i < n, < j < m, J\ k be the vector 

- n,m i m - oo,m 

space of fc dimensional vectors with entries in Y\ > and Ilm+i = Ilm+i' 

Utilizing the orthogonality relations (|3.1[) we obtain the following auxiliary re- 
sults. 

Lemma 3.1. Suppose $ € Jlfc • tf ^ satisfies the orthogonality relations, 

($, z l w j ) = 0, < i < n, < j < m, (3.9) 

then $ = T'lVm; where T is a k X (m + 1) matrix. If k = m + 1, T is upper 
triangular with positive diagonal entries, and if (<£>, $) = I m +\, then T = I m +i- 

Lemma 3.2. Suppose <£> G Jj. . //$ satisfies the orthogonality relations, 

(<£, z V) =0, < i < 7i, < j < 771, (3.10) 



f/ien $ = T$, ( 



where T is an k x (n + 1) matrix. If k = n + X, T is upper 



triangular with positive diagonal entries, and if (<£>, $) = I n +i, i/ien T = I, 



n+l- 



With the above we can make contact with the matrix orthogonal polynomials 
introduced in Section 2. This was observed by Delsarte et. al. [4] 



Lemma 3.3. Let Q n ,m be given by (|3.3p . Then 

x~ T m ( \\ m m— 1 

$n,m = L n (Z)[W , 10 



[1,10, ... , 7f m ] J m i?™(z) J, 



(3.11) 
(3.12) 
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<b n<m (z,w) 

$n-l,m{z,w) 



[w m ,w m -\...,l] T 



L 



Z J-m+l 
y n— 1 t 
1 J-m+l 



(3.13) 



Proof. If we substitute the equation 

® n , m = L n {z){w m ■■■ l] T = Y J Ln,z l [w m ■■■ if 

i 

into (13. 9p . where L n (z) is an (m + 1) x (m + 1) matrix polynomial of degree n, we 
find, for j = 0, . . . , n — 1, 



= ( $, 



1 



5> 

i=0 















1 




1 



Similarly, 



^ ' in,! 
n 

^i n! i£ m (z l , z 3 ) = C m {L n {z), z 1 ). 

i=l 

= C m (L n (z), L n (z)). 



This, coupled with (|2.8p and the fact that (|3.3p implies that L„ im is upper triangular 
with positive diagonal entries, gives p. lip . Equation (|3.12p follows from p. lip 
and (|2.36p while equation (|3 . 1 3|) follows from (|3.1ip and the definition of L. □ 

Analogous formulas for bivariate orthogonal polynomials in the reverse lexico- 
graphical ordering are obtained by interchanging the roles on n and m. 

The function M. given by (|2.3ip can be used to show that $ n . m satisfies a 
minimization condition. Define M : Y\. m +i ~ * Herm(m + 1) by 

M($) = ($t,$t)_($ + $t )) 

We find 

Lemma 3.4. The polynomial $ n ,m is the unique minimizer on Hm+i- 
Proof. Since $ £ Ilm+i can ^ e represented as 

n 

$(z, w) = [1, w, ...,w m ]$(z) = [l,w, ...,w m ] 

i=0 
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and from (l2~4l 



- w m- 




- w m- 


)■ 


" C NM {z l - ] ) ■ 


■ c NM (z^w- m y 


= C m (z\z i ) 


1 




_ 1 _ 




C NM {z^w m ) ■ 


£nm(z*~ j ) 





we find M(&) = M(&). The result now follows from (|3. 12|) and the fact that 
Rn,m R n( z ) minimizes M on Um+v 



□ 



We can now derive recurrence relations between the various polynomials. 



Theorem 3.5. Given {4>„ jm } and {& n ,m}, 0<n<N,0<m< M, the following 
recurrence formulas hold 



where 







(3.14) 






(3.15) 




r $ — $ i - JC <1> i 


(3.16) 






(3.17) 






(3.18) 




— Z" 1 cf> 4- (V 1 \ T< & T 
^ n,m J n,m ^ >V m i V x n,mJ n,ra— 1 5 


(3.19) 


E — 




(3.20) 


A n jyi — 


(z$ n _ 1 , m) # re r i JeM m + 1 - m + 1 , 


(3.21) 


^n,m — 


($n,m-l, $»-l,m) € M m ' n , 


(3.22) 


r — 


(*n >n ^l,*n,m) G M ro > m+1 , 


(3.23) 


JC 1 




(3.24) 


r 1 = 




(3.25) 






(3.26) 


i 1 = 

n , 7n 




(3.27) 



Similar formulas hold for <& n ,m and will be denoted by f 13.14|) -( 13.19l) . etc 



Proof. Equation (|3.14[) follows from Lemma [3~3l equations (|2.17[) . (|2.36p . and 
.Likewise (|3.15|) follows in an analogous manner from (|2.22|) . To prove (13.16 



2~37l) 



note 

that, because of the linear independence of the entries of $ njm , there is an m x 

a n— l,m 

(m + 1) matrix r„ )TO such that r„ )m $„ ]TTl - $„, m _i S rj m ■ Furthermore 

(r n , m $ n , m - $ n ,m-l,zV) =0, < i < n - 1 < j < m - 1. 
Thus Lemma 13.21 implies that 

r^m^n-m ^n,m— 1 — -^n.m^n — l,m- 

The remaining recurrence formulas follow in a similar manner. □ 
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Remark 3.6. As indicated in the proof, formula ()3 . 14(1 follows from the theory 
of matrix orthogonal polynomials and so allows us to compute in the n direction 
along a strip of size m + 1. This formula does not mix the polynomials in the 
two orderings. However, to increase m by one for polynomials constructed in the 
lexicographical ordering, the remaining relations show that orthogonal polynomials 
in the reverse lexicographical ordering must be used. 

Using the orthogonality relations from Lemma [37X1 Lemma [3~2l and (|3.ip we find 
the following relations. 

Proposition 3.7. The following relations hold between the coefficients in the equa- 
tions for $ and $. 



Also 



K — Aft T = ft 

f 1 - (J 1 ) T if 1 = (1C 1 ) T 

n.m v rt,m/ ' n.m \^n.m) s 



r r t = i - k, kl^ 

L n,m L n.m m '^n.m'^n.m, 
^n.mi^n.m) Im ^n.m (^n,m ) 

J rt _i_ ft p = j 

± n,m ± n,m ~ 1 n.m 1 n.m J m+1 

T 1 ( T 1 it i (pi itpi — j 

n,m\ n,m) ' \ n.m) n,m 1 m+l- 



(3.28) 
(3.29) 

(3.30) 
(3.31) 
(3.32) 
(3.33) 
(3.34) 



Remark 3.8. The matrix T rlym has a zero in the entries (i, j),i > j and has positive 
(i,i + 1) entries. Since r n , m rjj m = ^n^iU^UmT^^ where U m is the m x m + 1 
matrix given by 

U m = [0, I m ] , (3.35) 

we see that Y n ,mU m is the upper Cholesky factorization of the right hand side of 
(|3.31| . From this it is easy to obtain r n>m . The matrix m has zeroes in the 
entries > j with positive (i,i) entries. The matrix I n _ m has first row and 

column equal to zero except for a 1 in the (1,1) entry. 

The above recurrence formulas also give pointwise formulas for the recurrence 
coefficients. In order to obtain these formulas we define the m x m + 1 matrix f/i, 



U m = [I m , 0] , (3.36) 

and the (n + l)(m + 1) x (n + l)(m + 1) matrix P^'™ 1 which takes monomials in 
the lexicographical ordering to those in the reverse lexicographical ordering, i.e., 



P% m [z n w m ,z n w m - X , 



l} 1 = [w m z n lW m z n ~\...,l} 
Analogous equations hold for the n x [n + 1) matrices U n and 
Proposition 3.9. Let 



$ n ,m(z,w) = $™( Z ) 



and $ n , m (z,w) = $™ (w) 



(3.37) 



(3.38) 
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where 

$™(z) = $™ n Z n + $™ n _ 1 Z n - 1 + ..., 

= &L m ™ m + *™, m -i^"~ 1 + • • • , (3-39) 
then the following relations hold, 

r„, m = $™; 1 Z7 ro ($™„)- 1 , (3.40) 

r 1 = f$ m y 1 f3 4il 

n,m n,n u m\ n,n/ > V / 

^n,m — ^n,m^n,m^n,m^ (3.42) 

4,m - KJr^W.O, . . . ,0]C-^P;T T [^ + i,0, ■ • ■ ,0] r (^J-\ (3.44) 

^n,m = ( < J ) ™n T )~ 1 [0, ■ ■ ■ i 0, J m+ i]C Jl? 1 Tl F ) !] ,mT [7„ + i, 0, . . . ,0] T ($m,m)~ X i (3.45) 

where F = <&" r/ T ('<I>™ _ M _1 and F 1 = <l" (7/ 1, l T ('<E> ,l ~ 1 "l" 1 

Proof. Equation (]3 .41[) follows by equating the coefficients of z" in (|3 . 1 7[) on the left. 
The same argument gives (|3.41|) . To show (|3.42|) multiply (|3.18|) on the left by T n>m 
then subtract the resulting equation from (|3. 16|) . Now equating the coefficients of 
w m gives the result. Equation (|3.43[) follows by taking the transpose of the reverse 
of equation (|3.17[) then multiplying (|3.19p on the left by f x n m and subtracting the 
resulting equations. Equating powers of w m then gives the result. Equation (|3.44|) 
follows by equating the highest powers of w in equation (|3.18p and equation (|3.45|) 
follows in a similar manner from (|3 . 19[1 and the fact that C n ^ m is a doubly Tocplitz 
matrix. □ 

Remark 3.10. From (|3.11j) and Lemma 12.51 we see that ($™„)^ is the lower 
Cholesky factor of [/ m +i, 0, . . . , 0]C„ im [/ m +i, 0, . . . , 0] T and a similar relation holds 
between $^ m and C„ !m . Thus equations (|3.44|) and (|3.45|1 give the relation between 
In,m and l\ m and the Fourier coefficients of Cn.m- These coupled with equations 
(|3.42jl and (|3.43ll relate the Fourier coefficients of Cn,m to JC n ^ m and /C* m . 

We now give relations between the coefficients in the recurrence formulas at one 
level in terms of those at previous levels. 

Lemma 3.11 (Relations for /C„ jrra ). For < n, m, 

^n,m-l^n,m = ^n,m-l (^4„_i, m )^ — ^-"ra,m — l,m (-^n — 1 ,m )^ ' (3.46) 

^n,m(r„_i. m )' = A nm _iK, n —\^ m — A nm _ 1 E n _ m ^iiC n _ lm . (3.47) 

Proof. To show equation (|3.46[) multiply (|3.22[) on the left by T\ im _ 1 then use 
(|3.17p with m reduced by one to obtain, 

(u><I>„, m -2,$n-l .ml • 

Eliminating $„_i, m using ( 13TT3|) then applying (|3T22"|) and (|3~M)l gives IpTlS]) . 
Equation (|3.47|) follows in a analogous manner. □ 

Lemma 3.12 (Relations for K,^ m ). For < n,m, 

^n,m.-lfc n , m = ^■n.m-l{A n _ lm ) — K, n , m -\Fj n - i, rn (A rl _ 1 , m ) , (3.48) 
^-n,m(I"— l,m) = ^n.m-l^n-l.m ~ -^n.m— l^n,m— l^n— l,m- (3.49) 
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Proof. To show (|3~48l) multiply (I3T24]) on the left by r„, m _i then use ([37T6]) to 
obtain 

Now use ( 13.14j) with n reduced by one and then equations (|3.24p and (|3.22|) yields 
(I3.48p . Equation (|3.49|) follows in a similar manner. □ 

Lemma 3.13 (Relations for B n>m ). For < n,m, 

Tn— l,m-^n,m -^-n,m— \^n, m (-^n— l,m) ^n,7n— l^n — l,m' (3.50) 
-E , ri,m(r r i-l,m) = ^n-l,m(^ n ,m) ^n,m~ 1 + ^1- lm E nm _ i . (3.51) 

Proof. To establish (l3~50l) multiply |3T20|) on the left by r n _i, m then use (|3~T6)) to 
obtain 

With the use of (]3 . 14[) to eliminate z$ n _i. m _i we find 



.TTi— 1 ' n— l,m/* 



n— l,m) ^~ P"n,rn— 1 ( ^ n _ \ Ta _ i , 5* 5 

The second inner product on the right hand side of the above equation evaluates 
to r^ —1 m while the first may be evaluated using (|3.19[) followed by (|3.22p to give 
the claimed equation. To obtain (|3.5ip multiply (|3.20[) on the right by (T} 1 _ 1 m ) T , 
then use (|3.17p to get 

^n,m(r n -Lm) = ( z &n-l,m> ^n-l,m-l)- 

Using (|3.14p to eliminate $^-1 m-i yields 

En,m{Tn- 1 . m ) T = { z $n-l,m, $ n,m-l) ^n,m-l + (*n-l,m, ^n-l.m-l)^-! ■ 

Equation (|3~23l) can be used to evaluate the second inner product on the right hand 
side of the above equation while (|3.19p . ( 13.22p and (|3.28p can be used to obtain 
the first inner product. □ 

Lemma 3.14 (Relation for r* m ). For < n, m, 

^n,rnX\i,rn = ^n,m—l^n,m(I n ,m-l) ^n,m— l^Vs.m-l (3.52) 

Proof. To show (|332|l multiply (|3~23j) on the left by r+ m and use f3~TC)l to find, 

r n,m r n,m = (^n.m-l, $n,m-l) ~ (w$ n ,m-l) $n-l,rn)£l, m - (3.53) 

Eliminating u;$„ jm _i in the second term on the right hand side of the above equa- 
tion using (|3.17p then applying f 13.15p gives the third term on the right hand side 
of (|3.52p . In the first term on the right hand side of the above equation substitute 
the reverse transpose of Q3.19P to find 

{W$ n ,m-l,®n,m-l) = {w$n,m-l, ®n,m-l)( I n,m-l) T + T l,m-l T l,m-l> 

where (|3.23p has been used to obtain the second term on the right hand side of the 
above equation. The result may now be obtained by applying (|3.18p for $„ )m _ i 
and then using ( 13~20)1 □ 
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Lemma 3.15 (Relations for I n . m and 7„ m ). 

In,m^\i,m ^n,m^n,m] (3.54) 

^n,m = ~A n ^ m E n m A n ^ m I n ^ m + A n m I n _ l m T n ^ m , < n (3.55) 

Proof. Equation (|3.54[) follows by multiplying (|3.26p on the right by rj, m then 
using C I5TB . (l3~m and (l3~23l . In (l3~771) use f l3~T5D and (j3"25)l to find 

^n,m = ( ^ n,m> *&n,m)In,m + ( *£" 71, m > ^n— l,m)In,m • 

The first inner product on the right hand side may be evaluated using ()3.15p . To 
evaluate the second inner product eliminate $J lm using the reverse transpose of 
equation (|3. 15[) then use (|3.27|) to obtain the claimed equation. □ 



4. CHRISTOFFEL-DARBOUX FORMULAS 

The Christoffel-Darboux formula plays an important role in the theory of one 
variable scalar and matrix orthogonal polynomials. Using the connection between 
two variable orthogonal polynomials and matrix orthogonal polynomials we de- 
rive two variable analogs of the Christoffel-Darboux formula. These will play an 
important role in the theory of two variable stable polynomials discussed later. 

Lemma 4.1. Given {& n ,m} and {3> n ,m}, 

$„, m (z,<t/;) $^(21,^1) - z 1 z^ T nm {z,w)^\ m {z- l ,w 1 ) T (4.1a) 
= (1 - ziz)<5> n<m (z, w) T $r n ^ m {zi, wi) T 

+ $n-l,m{z, W)^l_ hm (zi, Wx) - w)*n-l >m (zi, Wif (4.1b) 

= (1 - z 1 z)$ n , m (z, w) T $l l>m (zi, wi) T 

+ ^n,m-l(z, W)% n^xiz^Wxf - ZiZ^ n>m _i(z, Vj) T $1 m _ j (z X , W 1 f . (4.1c) 

Proof. The equality (14.1al) ~ (|4.1bj) follows by subtracting (|2.24[) with n reduced by 
one from the original equation then using Lemma 13.31 The equality (|4.1a|) ~( 14.1c|) 
can be obtained in the following manner. Let 

Z n ,m{z, w) = [1, w, . . . , w m ][l m+1 ,zl m+1 , . . . , z n I m+1 ], 

and Z n , m (z, w) be given by a similar formula with the roles of z and w, and n and 
m interchanged. Then from Lemma [2751 (|2.24[) . and (|3.11|) we find 

$n,m{z, w) < S>i^ m {z ll W 1 ) - Z x Z^ m (z, w)& n m (Zj , W 1 ) T 
1 — Z\Z 

= Z n . m {z, w)C~^ m Z n , m {z\,W\)^ = Z ntTn (z, w)C~J m Z nim (z ll wi) f 

= $n,m( z i w )®n,m( z l> w l) T + Z n : ,n-1 (z, w)C~^ n _ 1 Z n , m _i {z X , Wi) f . 

Switching back to the lexicographical ordering in the second term in the last equa- 
tion then using Lemma 1^751 yields the result. □ 

As an immediate application of the above lemma we obtain, 
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Theorem 4.2 (Christoffel-Darboux formula). Given {&n,m} and {§n,m}, 

$n,m(z, W) $^ m (zi, Wl) - Z lZ ^l m (z, w)^j m {z 1 ,W 1 ) T 
1 — Z\Z 

n 

m 
3=0 

In the first line of the above equation the terms Z\z may be replaced by w\w if 
we switch to $ n ,m- 

An interesting variant of equation (|4. lc|) is, 

Lemma 4.3. 

$n,m(z, w) T $i, m (zi,W 1 ) T - C, m -l(^ w )^ n ,m-x( Z ^ W lf 

= $„, m (z,w) T $^ m (zi,wi) T -l>^_ lim (z, W )l>L l!m (zi, W i) T . (4.2) 

Proof. Equating the sums in the above Theorem yields 

m— 1 

ri-1 

= § n , m (*, w) T $l m (z uWl ) T - w ) $ J,mOl>™l) T - ( 4 - 4 ) 

Switching to the lexicographical ordering in the sum on the left hand side of the 
above equation and reverse lexicographical ordering in the sum on the right hand 
side, extracting the highest terms, then using the Christoffel-Darboux formula to 
eliminate the remaining sums gives the result. □ 

Remark 4.4. The above equations can be derived from the recurrence formulas in 
the previous sections. However, the derivation of equation (|4.1c[) is rather tedious. 

5. Algorithm 

In this section we use the relations developed earlier to provide an algorithm 
that allows us to compute the coefficients in the recurrence formula at higher levels 
in terms of those at lower levels plus some indeterminates that are equivalent to 
the moments. This will allow us to construct positive definite doubly Toeplitz 
matrices. As a byproduct we construct the orthogonal polynomials associated with 
these matrices. More precisely at each level we use the new indeterminates and the 
coefficients on the levels (n,m — I) and (n — I, to) to construct K, n ^ m and TC\ m . 
With this we can construct the other coefficients needed to proceed to the next 
level. The E n . m are closely related to the matrix recurrence coefficients needed to 
compute C n>m . Furthermore and & n ,m can also be computed. In order to 

construct the above matrices we will have need of the m x (m + 1) matrices U m and 
given by (|3.35[) and (|3.36p respectively, and the vector e™ € R m which is the 
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vector with one in the first entry and zeros everywhere else. From the definition of 
K} n m we see that 



fC„ 



(tfl$n,m-l, *n-l,m) 



~ jjn—X 
n,m— 1 



^n,m-l 



J n-1 
n— l,m 



Z ™ n-l,m 
j rnf n\T , r> 



(5.1) 



where R n ,m is an m x n matrix containing moments c$ j,{|i| < n, |j| < m} \ 
{(=bn,=bm)}. Likewise with the help of (|3.38[) and its tilde counterpart we find 



X,m) 





"0" 




"0" 




< — n,m^ n ,ro— 1 




v m,n— 1 




}t + i? 




1 




1 





(5.2) 



where R n ,rn contains only moments from lower levels. 

We proceed as follows, at level (0, 0) we have the parameter uo.o > which 
corresponds to cq o- The polynomials <I>o o and $o o are chosen as — . From (13.261) 

V u 0,0 

and (|3.27p we see that 7q,o = 1 — o- At l eve l (h 0) there is one new parameter w i0 
which can be taken to correspond with the one dimensional recurrence coefficient 
i.e. Uifi = <Xi = Ei t o corresponding to the (i, 0) level and must be less than one in 
magnitude. From (|3 . 30[) and the normalization chosen for the polynomials A^q = 

\Jl — l-E^ol 2 - This allows us to compute $i,o, and $,,o- The sizes of the matrices 
given in (j3~22f , (|3~23l) , (j3~24l) . and (f3T25|) show that 

fci,a = ^-i,o = r^o = K. i0 = K. i0 = T i0 = 0, 



where p.28[) , and (|3.29p have also been used. Furthermore (|3.33|) and Remark 
imply that — (e^ +1 ) T = lj where (13.28)) has been used. Equation f 13.3ip 
implies that I^o = Ui while (13.551) and (|3.29|) allow us to compute, 



E 



t 

i,0 



Ai.aE, 



■3,0 



(5.3) 



$i t o can now be computed from equation ( 13.18p . 

At level (0, j) there is one new parameter uo.j which as above can be taken to 



correspond with the one dimensional recurrence coefficient i.e. uq 



a 



Eq, 



corresponding to the (0,j) level and must be less than one in magnitude. The 
analysis for the (i, 0) level can be carried over with the roles of the lexicographical 
and reverse lexicographical orderings interchanged. Thus from ( 13.30[) and the 



normalization chosen for the polynomials Aqj 
compute $o,i- Again 



1 — l-Efljl 2 which allows us to 



^oj — Ko,j — foj — K,qj — K, j — f oj — 0. 
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Likewise Ioj — {e{ +1 ) T = Iqj and Tqj = Uj. Equations ( 13.55j) and (j3.29|) allow 
us to compute Jg j as above with i and j interchanged as well as the orderings. 
Equation (f3~18l) now allows us to compute $o j ■ 

At level (n, to) with n, m > there are two new parameters u n m and u_„ im since 
M -n,-m = u ntm and U n> - m = u-n,m- These along with the coefficients on (n — 1, m) 
and (n, m — 1) level will be used to compute K. n ,m and ^nm- This will be sufficient 
to compute the remaining coefficients on level (n., to). We begin with, 

Computation o//C„ !TO . If n = l,m = 1 then (|3.22|1 shows that ZC^i is a scalar 
which we choose as ui,-i- If m > 1 we see from (|3.41|) and (|3.36|) that 

pi 1 $™- 1 e m = 

n,m—l n,n J '< 

where e™ is the m-dimensional vector with zeros in all its entries except the last 
which is one. Since < &™~ 1 = A~ m _ 1 . . . Af m _ 1 $ 7l ~ is an upper triangular invert- 
ible matrix we find 

pl fam-lrrjjl \Ttt1 i m/ m\T\ _ pi S>m-l(Tjl \Ttj1 

and from (pT4Tj) ri m „ 1 $™- 1 (C/^_ 1 ) T = $™~ 2 . Thus equation (f3~46|) can be writ- 
ten as 

^-i(^ 1 )" 1 /c„. m ((^)t)- 1 

= (^n^n ) (£«,ra-l(^n-l,m)^ ~ £n,m-l^l-l,rn(^ri-l,m)^)((^m,m)^) 

= (^n,n ) (^Wi,m— 1 — ^■n,m—l^n—l,m) ((^m— l,m— l)^) 
— Hn,m—1- 

In the last equality we have used the fact that In-i,™!™"^ = ^^li m _i- Likewise, 

(^ 1 )- 1 /c„, m ((^)- 1 ) t (^-i) T 

= (*n-l,7i-l) H^n-l.m - -En,m-l^-l, m )((^mTm) ) 1 
— Hn— l,m* 

If 

(e™) T (^ 1 )- 1 /C„, m (fe t i)t)- 1 e^ = S n ,_ m , (5.4) 
then /C n ,m can be solved for as, 

(5.5) 

A necessary condition in order to be able to continue is that ||/C„ im || < 1. 

Computation ofT n ^ rn . Since /C n , TO is presumed to be a contraction Remark 13.81 
shows that L„ !m and L„ !m may be computed from the upper Cholesky factor of 

I - fcn,m%-n,m and 1 ~ ^n,m^-n,m respectively. 

Computation of IC^ m . In /C„ m we see from (j5.1[) that the only new entry is 
(£71,771)1,1 • If Ji = 1, to = 1 set /C1.1 = 1*1,1. If m > 1 we will show that all the rows 
except the first can be obtained from equation (|3.48|) . The structure of r„ iTra _i 
implies that r„. m _ie™ = so that 

r n ,m-i = r„ i „ i _i(L r ^_iJ7„ l _i +e™(e™) T ) = r nim _i?7^_ 1 t7 m _i. 
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But T niTn -iU^ l _ 1 is an invcrtiblc matrix, which allows us to rewrite (|3.48|) as follows 

U m — lK, nm = (rV^m— \U m _ i) {K"n,m— lC^n— l,m) — ^"n,m— \E n _ l,m(-^n-l,m) )• 

This gives all the entries in m except the first row. 
Similarly, if n > 1 we can write 

r„_i, m = r„_i >m (?7j_ 1 J7 n _i + e"(e™) T ) = r n _i )m u^ , _ 1 c/„_i, 

i.e., T^-i m = Un-i(Un-lFn_i m ), and equation (|3.49|) can be rewritten as 

^n,rrJ^n—X = iA«n,m—l^"n— l,m ~ -^71,711— l^n,m— \K"n— l,m){U n — l^n— l,m) (5-6) 

Thus the m x (n — 1) matrix /C^ m U^_ 1 , which is obtained from K} n m by deleting 
the first column, is known from the previous levels. This allows to compute all 
entries in the first row of K} n m except (/C* m )i,i and we put, 

(^Ln)l,l = u n ,m (5.7) 

A necessary condition on the parameters in order to be able to continue is that 
ll^-n mil < 1 which implies that |u njm | < 1 
Computation of E n ^ m 

We begin by taking the transpose of (|3 . 5 1 1) using the fact that E n ^ m is symmetric 
then multiplying on the left by the matrix e™ +1 (e™) T - Now multiply (I3.50|) by 
and add the resulting equations. If r n _i jm the (m + 1) x (to + 1) matrix obtained 
by stacking the first row of r„_j m on r„_i iTO we find 

fn-l,m^ri,m — ^7iG^-n,m— l^n,m (-^71— I.tti) -^7i,m— l-^n— l,7n) 

n,m—X^'rl,m^n—l,m ^n,m-lf n-l,m)- (5-8) 

From the structure of T 1 and T we see that r„_i, m is an upper triangular matrix 
with positive diagonal entries hence invertible. Thus E nm can be computed from 

the above equation. If H-EV^mH < 1 then E nyTn may be computed from f 13.50p and 
( 13.5ip . We may also compute A n ^ m , A n ^ m and the polynomials and $> n>m - 

While the condition that E n ^ m be a contraction is necessary and sufficient to be 
able to continue it is not optimal in the sense that it does not take into account the 
redundancy inherent in the equations giving E n ^ m . This will be taken into account 
in the computation of Y\ m . 

Computation o/r* m . As above we see that (|3.52[) gives 

^n,mU m (5-9) 

= {In,m—lE n ,m(I n ,m—l) + ^ n>m —lX n ,m—l 

^n,m^n—l,m^n—l,m^n—ltm^n,m){Um^ti,m) > 

which allows the computation of all the entries of m except the (1, 1) entry. Since 
(ef) T / n , m _i = (ef ) T , (ef ) T 4 im _ 1 = and likewise with / n , m _! and T\ i m _ 1 
replaced by I n ^ m and m -\ respectively we find with the help of (|5.8p . 

( e T) T ^n,mUm = ( e ™) T ^n,m + ( e lT) T ^-ri,m^ri,mJ (5.10) 
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where 

Hn,m = ^n,m-l((4,m-l)^n,™^n-l,ra + (Ifi )TO -l) -Er— l,m) (5-11) 

and 

Hn,m = ^n-l,m e l ( e l + ) (f^m-l) (-^n,m-l) (^mrj,,m) (5-12) 
+ ^n-l,m^n-l,m^n—l,m^n tm (U m J'n tm ) 

Thus the first entry T 1 can be computed using the first row of equation (|5. 10|) 
and (|3.32p which gives 

\Kmh,i)\ 2 = l-(e?) T Hl m (en (5-13) 

where 

rj-3 _ / tt2 i v-l ttI i IT 1 IT 1 "it I fr"l fK - ! \t ('c: i 4*1 

n,m V n,m ' '^n^m n,m)\ n,m ' ^n,irt n 7 mJ ' '^n,m\ n,m/ ' \ / 

Computation of the remaining coefficients. Using the arguments above we see 
that the relevant part of / n . m may be computed from (|3.54[) and l\ m may be 
computed from (|3.55|) . The matrix m can be computed in the same manner as 

r 1 



n,m " 



6. Construction of a positive linear functional 

The above algorithm allows us to find a linear functional given the coefficients 
in the recurrence formulas. More precisely, 

Theorem 6.1. Given parameters Uij £ C < i < n, \j\ < m,u-i_j — we 
construct 

• scalars Ei t o, i = 1, . . . ,n, and Eqj, j = 1, . . . , m; 

• i X j matrices /Qj i = 1, . . . , n, j = 1, . . . , m; 



• i x j numbers (e{) T Hf e{ i = 1, . . . , n, j = 1, . . . , m; 



If 

u ,o > 0; \Ei, \ < 1, |i ,il < 1, < 1, and <;'//;','! < 1 (6.1) 

then there exists a positive linear functional C on FJ"' m SU ch that 

£($i,m$lj = k, 3 I m +i and C(^ n ^l d ) = 5 id I n+x . (6.2) 
The conditions (|6.1|) are also necessary. 

Proof. We construct the linear functional by induction. First, if n — m — we set 
£(1) = uo,o and $ ,o = $o,o 



/uo,o 

and thus £($ ,o$o,o) = £($o,o*o,o) = L 



If m = 0, we construct Ai,o = y 1 — \Ei,o\ 2 where E^.o = Uj,o- The polynomials 
$1,0 i — 0, ... n are now computed using (|3. 14|) and then we define 

£(*i,0*}, ) 

This gives a well defined positive linear functional on z J for |j| < n. 
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Likewise, if n = 0, we construct $o,fc using ( 13.14j) and define 

C(<f>o,M,j) = 

which gives the linear functional on w J for \j\ < to. Thus formula (16. 2p will hold if 
to = or ?i = 0. 

Assume now that the functional C is well defined and positive for all levels 
0<i<n — 1, 0<j<m and < i < n, < j < m — 1 before (n, to). To ease 
notation we will use the bracket given in equation (|3.8p with £n,m replaced by C. 
We first extend C so that 



-l,$n-l, 



/Cr, 



(6.3) 



To check that the above equation is consistent with how C is defined on the previous 
levels, note that from (|3.46l) 



,m—l^n,m— 1j ^n— l,m) — m _ l^n,mi 



(6.4) 



which follows from the construction of /C„ im and the definition of £ on the previous 
levels (see Lemma l3.1ip . Similarly, using the second defining relation of /C„ iiri (i.e., 
the last row of ([3747])) we see that 



^-n,m(r n — 1 



It 



(6.5) 



Equations (|6.4p and (|6.5p show that most of (|6.3p is automatically true. We now 
define C(z n w~ m ) so that (|5.4|) holds which completes (|6.3p . 

Using an analogous argument we can use the construction of fC„ m to extend the 
functional to z n w m so that 



(6.6) 



This completes the extension of £. What remains to show are the equations (|6.2p 
hold. This is accomplished by first constructing E n-m from H5.8P . The condition 
on (e™) T -ff^ m e ]™ an< i l|5.9p show that the first row of m may be computed and 
that we may choose 

( r n,m)l,l > 0. 

With the first row of r* and all of r„ !m (which is calculated from the Cholesky 
factorization of /C nim /C£ m ) $ n ,m may be constructed from (|3.16p and (|3.17p . Equa- 
tions (|6.3p and (|6.6p coupled with (|3.16p . (|3.17p and the orthogonality relations on 
the previous levels show that 



(r„ m $„ >m , $„_i lfc ) = 0, fc = 0, 1, 



and 



((er) T ri m $„, m! w fc 

Equations (|6.7p and (|3.17p show 



o = ((ef) T r* $„, ro , 



> = o, 



. ,TO, 



1,2,. 



(6.7) 



i) — (( e i , ) T rJ im ^ti,mi 
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- Z n-Y 

The fact that &n-im ^as an invertible coefficient multiplying 



has been 



f 

used to obtain the second equality in the above equation. The above implies 

($n,m, $n-l,k) =0, fc = 0, 1, . . . , TO, 

which in turn implies that 

($n,m, ®j,m) = 0, j = 0, 1, . . . , Tl - 1. 

To show that 

we note that equations (|3.16p , (|3.f 7p and (|3.52[) imply 

/fpm^r 1 a> r * \ — (p m \ T v 1 r+ 

and (|3.32jl implies 

((e™) T r,\ m <£y m , (e™) T r^ lTO <i> nm ) = (e™) T r^ m (r^ m )^(e™). 

Thus £ is a positive linear functional. The orthogonality relations for the polyno- 
mials $>ij now follow. □ 

Let C(T 2 ) denote the set of continuous functions on the bi-circle, above Theorem 
now allows, 

Theorem 6.2. Given parameters Ui : j S C with Ui —j = Q-jj- If equations (|6.ip 
hold for all < i,j then there exists a positive measure /i supported on the bi-circle 
such that for any f G C(T 2 ), 



T- 



Proof. From the hypotheses imposed above, Theorem 16 . 1 1 shows that C„ jm is pos- 
itive definite for all n and m so the result follows from Bochner's Theorem jT8j 
Section 1.4.3]. □ 

Remark 6.3. The above construction gives a criteria for the existence of a one 
step extension of the functional. That is, given moments so that there exists a 
positive linear functional on FJ" ' m U FT™'"* - ; an y se ^ 

{^n,mi n.m\ H—n. — m — ^n.m ^n, — m — n.ra 

that satisfies (|6.1j) can be used to extend the functional to FJ ra '"\ However it is not 
difficult to construct examples where no extension exists. See section 8.2 

7. Two Variable Stable polynomials and Fejer-Riesz factorization 

In this section we study the consequences of K, n ,m = 0. This will make a connec- 
tion with the results in [8] on stable polynomials and the Fejer-Riesz factorization 
theorem. We say that a polynomial p(z, w) is stable if p(z, w) ^ 0, \z\ < 1, \w\ < 1. 
A polynomial p is of degree (n, m) if 



i=a j=o 
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with k n ,m 7^ 0- Finally we say that the polynomial p n ,m of degree (n,m) has the 
spectral matching property (up to (n,m)) if 

C i zkw °) = 7^2 / i *f M2 m & z = e ie ,w = e^, 
(2iry J T 2 \p n , m ( z ' w )\ 



for |fc| < n, \j\ < to. 

Lemma 7.1. Si 

K-n, m = 0, then, 



Lemma 7.1. Suppose that C is a positive definite linear functional on rj"' m ano 



4> n,rn( z > w ) <P n,m( Z ^ w l) ~ K]rn( z > w ) ( f>n,m( Z U'W 1 ) 
= (1 - WWi)$ n , m -i(z, w)*L m _i(zi, Wi) 

+ (1 - ZZi)$ n - ljTn (z, w) T ^ n _ hm (z 1 , Wi). (7.1) 

Proof. If K n;m — then (|3.16|) shows that T rhm <& rh . m {z,w) — G> ntm -i(z,w). Thus 

$n, m (z, w)Tl m = w$ ntm -i(z,w). Also (|3.31|) implies that (r„, m )(i ji+ i) = 1, i = 
1..to with all other entries zero. Thus we find 



( l > n,m( z , w )'Pn,m( z l, w l) + $ n,m {z, Vj)^ tTn T n , m $ n<m (z 1 , Wl ) t 



= < frn,m( Z > W ) <f>™m( Z li W l) + ™l * ri,™-l (z, ») $ „, m -l (21, Wl ) + . (7.2) 

From (|4.1c|) we find 



<f>n,m( z ' w ) 4>™ m ( z l, Wl) ~ ZZl<f>™ m (z,w)ff)™ m (zi,Wi) 
= (1 - WWi)$n,m-l(z,w)$n, m -l(Zl,Wl) 
+ (1 - zz~i)$„ tm (z, w) T ^l i m (zi,Wi). 

Using ( 13.16P and the fact that ^ m (z,w) = <f>™ m gives the result. □ 
We now have 

Theorem 7.2. Suppose that £ is a positive definite linear functional on Y[ n ' m and 
K"n,m = then (f>™ m (z,w) is stable and, 

1 f p —ik8p—il<f> 

Conversely if ■n n _ m (z,w) is a polynomial of degree {n, to) suc/i that 7T„,m * s stable 

, I f e -ike e -U4 

Jr> k n ,m(e i9 ,e'*)|2' 
i/ien K. n ^ m = 0. 

Proof. If £ is positive definite and /C„ !m is equal to zero then Lemma l7.1l shows that 
(j)™ m (z,w) satisfies (I7.ll) . The first part of the result now follows from the proof of 
Theorem 2.3.1 in [5]. To show the second part let f(z,w) = 1/|V n , m {z, w)\ 2 , \z\ = 
1 = \w\ be the spectral density function associated with ~n n ,m- Then from equa- 
tion 2.1.5 in [8] and Lemma T3. 31 we find that 

$ n<m (z,w) = [V„, m (z, w), W$ n , m -i(z,w)]. 



1 f p—ikBp—iltf) 

C{e- ik6 e- il t>) = (-) 2 / ————d9d<P, \k\ < n, \l\ < to, 
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But this implies that $„ im (z, w) = [TT n ,m (z,w),<Jy m _i(z,u>) T ] T . Hence from (I3J61) 

fcn,rn — 0. EH 

This leads to the following alternative proof of the two- variable Fejer-Riesz The- 
orem in [8] . 

Theorem 7.3. Suppose that f(z, w) = J2k=-n YaL— m fki zku > 1 * s positive for \z\ = 
\w\ = 1. Then there exists a polynomial 

n m 

p(z,w) = ^2^2p u z k w l 

k=Q 1=0 

with p(z, w) ^ for \z\, \ w\ < 1, and f(z, w) = \p{z, w)\ 2 if and only if tC ntm = 0. 

Proof. For g E C(T 2 ) let C{g) = J T2 lp ffi*l ){2 d6d(f>. Then £ is a positive 
definite linear functional on T 2 . The necessary part of the above result now follows 
from Theorem 17.21 The sufficiency also follows from the above Theorem and the 
maximal entropy condition [T]. □ 

An alternative approach for finding a factorization as above may be done using 
the notion of intersecting zeros (see [9]). Also, the question of factorizing a non- 
negative trigonometric polynomial as a modulus square of an outer polynomial was 
addressed in [B], allowing for generalizations in the operator valued case. When 
such a factorization of the desired degree does not exist one can approximate the 
trigonometric polynomial with one that does have the desired factorization. This 
question was pursued in [12] . 

The vanishing of K, n ^ m has the following geometric interpretation. 

Lemma 7.4. Suppose C is positive definite on FJ ' then K, n , m = if and only if 
for $„ !m _i constructed as in (|3 . 3[) , 

($ n ,m-i,zV) =0 j = m,0<i<n-l. (7.3) 

Proof. The definition of <&„ !m _i shows that it is already orthogonal to z % uP , < 
i < n, < j < m — 1. The remaining orthogonality conditions show that $„ !m _i 
is orthogonal to all the monomials in $„_i. m . Thus the sufficiency part of the 
Theorem follows from (|3.22p . To see the necessary part note that from the definition 



lCn,m — ($n,m-l,$n-l,m 



(7.4) 



with m an invertible matrix. Thus (|7.3p follows. □ 



Unfortunately at this point we are unable to see what the condition JC n 
implies for Uij, \i\ < n, \j\ < m except for u n ,- m = which follows from 
equation (|5.4[) . We can however get a partial characterization for when a positive 
measure on the bi-circle can be written as the reciprocal of the magnitude square 
of a stable polynomial. We begin with the following auxiliary result. 

Lemma 7.5. If Ei.j = then the first column of K\ ■ is equal to zero, in partic- 
ular Uij — 0. If Eij, ICij and /Cj_ij are zero then so is Eij-i. Conversely if 
tCij, Ki-ij, Eij-i andui.j are zero then E^j — 0. In both cases K,}, = [0, A. 
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Likewise if Eij = then the first row of K}j is equal to zero. If Eij, fCij and 
fCij—i are zero then so is Ei j—i. Conversely if ICi_j, Ki^-\, Ei-ij and Uij are 
zero then Eij — 0. In both cases K.jj = [0, (JC\ j_ 1 ) T ] T ■ 

Proof. If Eij — then equation (|3.5ip , and Remark 13.81 show that the first column 
of JC] ■ is zero. If /Q-i^ is equal to zero then (|3.54|) shows that all the entries of 
are all zero except for a one in the first entry. Thus (|3.50|l and (|3.51j) imply that 
if E id = 0, fCi tj = 0, and JQ-ij = then ^■_if^_ 1J = and Eij^f^j = 0. 
Following the argument in the construction of E n ^ m we see that -Bij-i = 0. The 
above hypothesis on K-i-i.j shows that ^i-ij — Ui-i thus (|5.6p and the fact that 
the first column of K,\ ^ is zero gives K,\ ■ = [0, /C^_ x . The converse statement 
follows from equation (|5.8|) . The remaining statements follow in an analogous 
fashion using Proposition 13.71 . □ 

Lemma 7.6. Let /i be a positive measure on the bi-circle. Then /j, is purely abso- 
lutely continuous with respect to Lebesgue measure and 

dn(6, dp) = —^d6d<t>, 

\Pn,m I 

where p n ,m is a polynomial of degree (n,m) with *p n ,m{ z i w ) stable if and only if 

K-i.j = 0, Ei+ij = 0, and E n ,j + i = 0, i > n, j > m. 
Proof. Suppose that dfi = t- } ., 3 d6d(f> with < p n ,m stable, then the sequence 

{ipi,j(z,w)} ip it j(z,w) = z l ~ n w j ~ m p n>m (z,w), i>n,j> m, 

is a set of polynomials with degrees (i,j) respectively such that tp ij — *p n ,m are 
stable and have the spectral matching property. Thus Theorem 17.21 implies that 

JCi.j = for i > n, j > m. Since — ^i.ji i ^ n i 3 > 171 we see from 

equation 2.1.5 in [5] that = $jj for i > n , j > m. This implies that 

— Ij+i so that equation (|3 . 14[1 shows that -Ej+ij =0, i > n, j > m. Since 

4>\j + i — j, i > n, j > m the preceding argument shows that S-tj+i = 0, i > 
n i j > ™- This proves the necessary part. 

To prove sufficiency note that if ICi j = 0, i > n, j > m there exist polyno- 
mials ?pij of degree where is a stable polynomial which has the spec- 
tral matching property. In order to show that ip ij = %j) n ,m we note that since 
Ei+i.j — equation (13.141) implies that = ^i,j> * > 3>rn. Furthermore 

Enj+i — 0, j > m implies that $ n ,j'+i = $n,m- Since Vij = $y = 4>\j for 
i> n, j > ti the result follows. □ 

The conditions on /Ci.j, and -EVij given in Lemma 17.61 are not optimal since 
they are redundant. Some of this redundancy is removed in the next theorem. 

Theorem 7.7. Let /i be a positive measure on the bi-circle. Then [i is purely 
absolutely continuous with respect to Lebesgue measure and d[i = where 
Pn,m is a polynomial of degree (n, m) with *p n ,m stable if and only if 
a. K- n .j — 0, U„_i J+ i = and — 0, j > m, 

b- /C ijTO = 0, E i>m -i = and Ui, m = 0, i > n 
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c. uui j = 0, i > n, j > m. 

Remark 7.8. Equation (|5.4[) and Lemma 17.51 show that U- n ,j, u~i,m, Un—i,j+i 
and Ui t m-i are also equal to zero for j > m,i > n. 

Proof. If fi has the form indicated in the hypotheses then Lemma 17.61 says that 
JCi_j = 0, i > n, j > m which coupled with equation (|5.4|) implies that U-ij = 
0, i > n, j > to the remaining conditions on the coefficients follow from Lemma 

17.51 If the coefficients obey a.-c. then Lemma [7.51 shows that £^ + i iiri and E n ,j + i, 

are equal to zero for i > n and j > to. Since E nm +i — 0, E n+ i Jm = and by 
hypothesis u_„_i, m +i = 0, equation (|5.5[) shows that /C n +i, m +i = 0. With this 

Lemma 17.51 shows that E n+ i^ n+ \ = and £Vi+i,m+i = 0. The result now follows 
by induction. □ 

It is possible to modify slightly the hypotheses of Theorem l7.7l to obtain a state- 
ment just on the coefficients in the recurrence formulas. 



Theorem 7.9. Suppose Ui.j are given so that equations (|6.1j) are satisfied for < 
i < n — 171 an d so that a-c in Theorem \ 7. 7| hold. Then for f £ C(T 2 ), 

1 



'2tt' 



T 2 



where \x is absolutely continuous with respect to Lebesgue measure with d\i — m 

Here p n ,m *s a polynomial of degree (n, to) with *p n ,m stable. 

Proof. From Theorem 16.11 that there exists a positive definite linear functional on 
rj n,m with the above parameters and from Theorem 17.21 the functional has the 
representation 

i f -ikd -iltf> 

£(.-*«-*) . ( -r I W < n, VI < m . 

with Pn^m a polynomial of degree (n, to) with *p n ,m stable. The result now follows 
from Theorem 17.71 □ 



8. Examples 

We now give some examples that illustrate various aspects of the results pre- 
sented earlier. We begin with the case n = 1, to = 1 with uo,o = 1, ^Ci.i = = 
and K\ 1 — u\ y \. From Theorem 16.11 we see that we must chose |iio,i| < 1 and 
|ui,o| < 1. Since /Ci.i = the only remaining condition for £ to be a positive 
linear functional on H 1 ' 1 is for e{ T Hf x eJ < 1. From equation (|5.9p we see that 

\Ui = IifiE\.i {l\ ) T . The construction of J^o, I\q and equation (|5.8[) shows 
that e{ T Hf < 1 is given by 

a|"i,i| 2 + &0"i,i^o,i u i,O + ui.itto.iui.o) + c < 1, 

with a = 1 ~ l "?- 1 "%° |2 , 6 = Vi_±^ill and c = |u i| 2 . This simplifies to 
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where 



(1 - |tto,i«i,o| 2 )«i,i 

U\ 1 = — = — = + Un lUi o- 

1,1 vi-l^ilVi-Kol 2 



Thus from Theorem 16. II and Theorem 17.21 we see that with it .o = 1, 

C(e- ike e-^) = (^) 2 / u e "Hr!l 2 ^# 1*1 < 1 S lil < 1. 

where (f>i^i constructed using equation (|3 . 1 6[) and the top row of (|3.17p is a polyno- 
mial of degree (1,1) with 0i i stable if and only if |uo,i| < 1, \ui t o\ < 1, U-i,i = 
and \ui,i\ < 1- Furthermore if we set u^o, uqj, Uij equal to zero for i > 1, |j| > 1 
then Theorem 17.91 shows that the above representation for £ extends to all contin- 
uous functions on T 2 . 

We can also use the previous results to investigate contractive Toeplitz matrices. 
In this case we find 

" I CU] 
Ci I 



Cm = 



(8.1) 



where C_i = Cj is a 2 x 2 Toeplitz matrix. In this case wo,o = 1 an d u o,i = so 
that E ,i — 0, A Q _i — 1. Since /Ci^i = we find I\i = [0, ^Jl — | 2 ] . This 

plus the computation of T\ described in the construction of C yields 

/=(e 1 ) T ^ 1 3 1 (ff 1 3 il )t(e 1 ) (8.2) 
= (1 + d)|wi,i| 2 + d(ui,iti-i,i -I- ui,iu_i,i) + <i|u_i 7 i| 2 < 1, 

where 

Ko| 2 



(l-Kol^a-Kr.xl 2 )" 
By completing the square this can be simplified to 

\u\,i\ < 1 

where 

= (1 + d)Jl - |ui, | 2 ui,i + 

and 

l^l 2 



(l-ltt-uHv^KoP 

which puts constraints on it_i : i. Thus we find the conditions for £ to be a positive 
linear functional and hence G\ to be a contractive Toeplitz matrix are 1 < 1, 
\v,-i i| < 1 and i| < 1. These constraints may not be strong enough to allow £ 
to be extended. To sec this suppose n — 1, m = 2, u .2 = and u 10 = 0. It is not 
difficult to see then that E\ \ = diag(ui i, u_\ t \). With U\ o = the constraint on 
w^i above reduces to |iti,i| < 1. However 

/ ti-i,i \ 
/ (l- 

'1.2 



K = | (l-l^l.ll 2 ) 1 / 2 



\(i-|«_i,i| 2 )V2 y 

so we see that in order for i^i 2 to be a contraction I"- 1 - 1 ! < 1, which may not 
be satisfied. 
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